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*^j Abstract. We enumerate a certain class of monomino-domino 

coverings of square grids, which conform to the tatami restriction; 

no four tiles meet. 

Let T„ be the set of monomino-domino tatami coverings of the 
i— i nx n grid with the maximum number, n, of monominocs, oriented 

Q^ so that they have a monomino in each of the top left and top 

f) right corners. We give an algorithm for exhaustively generating 

the coverings in T„ with exactly v vertical dominoes in constant 
^3 amortized time, and an explicit formula for counting them. The 

polynomial that generates these counts has the factorisation 

where S n (z) = n_iLi(l + z% )' ano - Pn{z) is an irreducible polyno- 
f- ^ mial, at least for 1 < n < 200. We present some compelling proper- 

ty — ties and conjectures about P n (z). For example P n (l) = n2"(™ _2 ) _1 

^D for all n > 2, where v{ri) is the number of Is in the binary repre- 

^—; sentation of n and deg(P„(z)) = J2k=i Od(k), where Od(k) is the 

largest odd divisor of k. 

o 

CO 

1. Introduction 

S^ The counting of domino coverings, together with its extension to 

counting perfect matchings in (planar) graphs, is a classic area of enu- 
merative combinatorics and theoretical computer science. Less atten- 
tion has been paid, however, to problems where the local interactions of 
the dominoes are restricted in some fashion. Perhaps the most natural 
such restriction is the "tatami" condition, defined below. The tatami 
condition is quite restrictive: for example, the 10 x 13 grid cannot be 
covered with dominoes and also satisfy the tatami condition. In this 
paper we restrict our attention to square grids, and explore in some 
detail the enumeration of certain extremal configurations. 

Tatami mats are a traditional Japanese floor covering whose dimen- 
sions are approximately lm x lm or lm x 2m. In certain arrangements, 

l 
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no four tatami mats may meet. Such an arrangement has a preferable 
structure which is discussed in p] and [2]. 

A tatami covering is an arrangement of 1 x 1 monominoes, 1x2 
horizontal dominoes, and 2x1 vertical dominoes, in which no four 
tiles meet. The present discussion is about tatami coverings of the 
n x n grid with exactly n monominoes and v vertical dominoes. On 
the basis of some computer investigations, Don Knuth discovered that 
the generating polynomial for small tatami coverings of this type, with 
respect to the number of vertical dominoes they contain, is a product 
of cyclotomic polynomials and a mainly mysterious, irreducible poly- 
nomial (private communication, December 2010). Knuth's discovery 
and our own observations motivated Conjecture 4 in [1], which is pre- 
sented here as Equation (J2l). In this paper we generalize and prove 
Knuth's cyclotomic factors, and determine some important properties 
of the mysterious polynomial. 

We prove in [TJ that all of the n x n coverings with n monominoes can 
be rotated so that monominoes appear in each of the top two corners of 
the grid, so we let T n be the set of these. Let H(n, k) be the number of 
coverings in T n with exactly k horizontal dominoes, and let V(n, k) be 
the number with exactly k vertical dominoes. Let S n (z) = Yl^ii^+z 1 ). 
We prove that the polynomial 

(1) W n (z) : = 2 J2 S n ^ 2 (z)S^ 1 (z)z n - t ' 1 + (S^j (z)Y , 

is equal to J2k>o-H( n ,k)z k for odd n and is equal to J2 k>0 V(n, k)z k 
for even n. Knuth's observation generalizes to 

(2) W n (z)=P n (z)l[S^^z), 

where P n (z) is the "mysterious" polynomial. We prove Equation pi) 
in Theorem HI 

The remaining factors of Equation (J2J) are of the form Sk(z), where 
A; is a binary right shift of n — 2, and the complete factorisation of these 
is known in general. The ith cyclotomic polynomial, §i(z), is defined as 
ELef^ -2 — u )i wnere ^ is the set of ith primitive roots of unity. Lemma 
5 in pQ shows that Sk(z) is a certain product of cyclotomic polynomials, 
which are known to be irreducible, and thus VH n (z) can apparently be 
factored completely as 

(3) W n (z) = P n (z) 11 $ 2j (z) L^J . 
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We have verified the irreducibility of P n {z) for 1 < n < 200 (the de- 
gree of Pigg(z) is 13022 and its largest coefficient has 55 digits), and 
thus we hope that Equation (|3| is the complete factorisation of \H n (z) 
for all n > 2. 

The class P n {z) of polynomials has some compelling properties, some 
of which are theory, others empirical. For example, we observe in 
Conjecture [2] that the alternating sums of P n (z) are the coefficients of 
the ordinary generating function 

Ef,(-»-' = ft + *)(l-2f)_ 

for 1 < n < 200. If the conjecture is true, then P 2 (n+i)( — 1) = ( ")• 
Furthermore P n (— 1) is equal to the sum of the absolute values of the 
coefficients of P n (z), only for n > 20. This second fact is surprising, 
considering the way P n (z) is derived - why n > 20? 

The complex roots of P n (z) appear to cluster neatly around the unit 
circle, and form convergent sequences as n — > 00. They are plotted 
in Fig. [7] for odd n; for even n, the plot has a similar look. 

Theoretical progress on P n (z) comprises Theorem [4] and Theorem [5] 
The former states that deg(P n (z)) = Y^kZi Od(k), where Od{n) is the 
largest odd divisor of n. We prove in Theorem [5] that for all n > 2, 
the sum of the coefficients of P n (z) is equal to n2 v ^ n ~ 2 ^~ 1 , where v{n) 
is the number of 1-bits in the binary representation of n. 

Our technique for finding \H n (z) employs an operation which pre- 
serves the tatami condition, called the diagonal flip, defined in [2]. The 
added observation that a diagonal flip changes the orientation of some 
dominoes, enables us to further exploit it. The crux of the argument 
uses the partition of T n , from Theorem 2 of [2j, which reveals diagonal 
flips each with 1,2, ... ,k dominoes, respectively, that can be flipped 
independently. We use this to express \H n (z) in terms of Sk(z), the 
generating polynomial for the number of subsets of {1,2, ... ,k} whose 
elements sum to i. 

The formula for the fth coefficient of \H n (z) translates into an algo- 
rithm for generating all possible nxn tatami coverings with v vertical 
(or horizontal) dominoes, given that we have one for generating all k- 
sum subsets of the n set, for k,n > 0. We employ an algorithm from 
[7] to generate our coverings in constant amortized time. 

1.1. Overview. The structure of square tatami coverings with the 
maximum number of monominoes is summarized in Section [2] — see 
reference pQ for a complete proof. In Section [3] we describe a represen- 
tation for the coverings in T n using strings over a ternary alphabet. 
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Each symbol represents a monomino which can either be flipped in 
exactly one of two diagonals, or unflipped. This representation is es- 
sentially the same as the one in j2], but our revision of the indices 
greatly simplifies the counting of vertical and horizontal dominoes. 

Our main result is Theorem [IJ the generating polynomial \H n (z), 
but the technical part is proving the formula for its coefficients, which 
is presented in Lemma [I] We discuss the coefficients of Xff n (z) and 
derive the self-reciprocal generating polynomial required to partition 
the coverings in T n and their four distinct rotations according to their 
numbers of vertical dominos. 

In Section [5] we apply the proof of Lemma [l] to generate T n in con- 
stant amortized time by adapting an algorithm given in [7j. 

In Section M we depart from the geometrical interpretation of VH n (z) 
and prove Theorem [3j the factorisation in Equation pi), and the re- 
mainder of the section focuses on properties of the factor P n (z), in- 
cluding Theorem [4] and Theorem |5j introduced above. 



2. Structure of a tatami covering 

An edge refers to the edge of a tile, while a boundary refers to the 
outer boundary of a covering. 

The structure of coverings in T n is characterized in Lemma 3 and 
Corollary 2 of pp. Corollary 2 states that an n x n covering with n 
monominoes has monominoes in exactly two of its corners, which must 
share a boundary. This ensures that each covering of T n has distinct 
rotations through 0, tt/2, tt, and 37r/2 radians, and that rotations of 
distinct coverings in T n are distinct from each other. 

If T e T n , then a diagonal, D, of T is a contiguous sequence of 
like- aligned dominoes whose centers lie on a line with slope 1 or — 1. 
The sequence must begin with a domino with its long edge on the 
boundary; the final domino will share an edge with a monomino which 



is also considered to be part of the diagonal. In Fig. 1(a) there are 12 
diagonals, only one of which contains horizontal dominoes. 

A diagonal flip of D consists of removing it from T, reflecting hori- 
zontally, rotating by | radians, and placing it back onto the grid squares 
that were vacated. 

There are three things to note about the diagonal flip: 

• a flipped diagonal is a diagonal; 

• the operation preserves the tatami restriction; and, 

• it changes the orientation of the dominoes that it contains, and 
maps the monomino to the other extreme of the diagonal. 



x n TATAMI COVERINGS WITH v VERTICAL DOMINOES 



The running bond, or simply bond, is a rotation of the basic brick 
laying pattern, in which all dominoes have the same orientation. The 
restriction that coverings in T n have monominoes in their upper cor- 
ners, implies that exactly one bond pattern is possible for each n. When 
n is even, the bond consists of horizontal dominoes, with monominoes 
along the left and right boundaries, and when n is odd, the bond con- 
sists of vertical dominoes, with monominoes along the top and bottom 



boundaries (see Fig. 1(b)). 

Lemma 3 in pQ shows that every covering in T n can be produced 
from the running bond, via a finite sequence of diagonal flips in which 
each monomino is moved at most once, and the top corner monominoes 
are not moved at all. In addition, any such sequence of flips results in 
an element of T n (see Fig. 1(c) \ . 



(a) 



(b) 




Figure 1. Examples of coverings in T n . (a) A diagonal 
flip in a 9 x 9 covering, (b) The horizontal running bond 
for n = 10. (c) A sequence of diagonal flips results in a 
covering in T 10 . Flipped monominoes are coloured red. 

From this perspective, we can look upon the original position of a 
monomino as its place in the running bond, and then describe it as 
flipped in a given direction if it has moved from its original position in 
T; otherwise it is unflipped. 



3. Representing coverings as a ternary string 

We describe a ternary string representation for n x n coverings with 
n monominoes. Recall that each monomino, besides the two corner 
monominoes, is in exactly two diagonals in the running bond, and in 
a given covering a monomino is flipped in one of these diagonals, or it 
is unflipped. A ternary symbol for each monomino indicates which of 
the three possible states it assumes. Each covering is described by a 
unique string of these ternary symbols, represented in the same order 



as the following indexed labelling (see caption at Fig. 2(c)). 
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Monominoes and their diagonals are labelled as shown in Fig. [2j such 
that the index, i, of a monomino is equal to the length of one of its 
diagonals, and n — i — 1 is the length of the other. This relationship 
between diagonal length and index is helpful in Lemma [1} 

The ternary string representing the 10 x 10 covering in Fig. 2(c) is 
s = (0, 1, —1,0,0, 1, —1), where s, = 1 if the i th monomino is flipped 
upward, Sj = — 1 if it is flipped downward, and Sj = if it is unflipped. 
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Figure 2. Labelling for T„. (a) For odd n, monominoes 
are labelled £j and bi. The distances from tj and bi to the 
left boundary are both i. (b) For even n, monominoes 
are labelled U and r^. The distances from U to the bottom 
boundary, and from rj to the top boundary, are both i. 
(c) The covering, (0,1,-1,0,0,1,-1). 

We use li,lj,rj,rf,t^,tj~, b^ , b^~ to denote the diagonals that the 
monominoes li, r iy t iy bi can be flipped on. Naturally, li and rj can only 
be (diagonally) flipped up or down, whilst U and bi can only be flipped 
left or right. 

Let d n (a) be the number of dominoes in the diagonal a, also called 
the length or size of the diagonal. It is a function of the index and 
direction of a: 



dJa) 



i, iiae{lir},tt,bf}; 
n-i-1, iiae{lj,rf,t^,b^}. 



Flipped diagonals which intersect are called conflicting, and can oc- 
cur as one of two types (see Fig. [3]). 



x n TATAMI COVERINGS WITH v VERTICAL DOMINOES 



- -L-i- 


1 1 







(a) 



(b) 



Figure 3. Example of, (a), Type 1 conflict, and, (b), 
Type 2 conflict. 
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Table 1. Conditions for Type 2 conflicts. 



Type 1: A pair of diagonals with monominoes originating on the 
same boundary are flipped toward one another (e.g. (t^,t^~) 
for some i < j). 

Type 2: A pair of diagonals with monominoes originating on op- 
posite boundaries are flipped in the same direction (e.g. (/J, rj)) 
and their combined length is at least n (see Table fl]). 

Lastly, if a is a diagonal containing a given monomino, let a be the 
monomino's other diagonal. 

3.1. A partition of T n . Let T n (a) C T n , where a is a diagonal such 
that d n (a) > d n (a), be defined as the collection of coverings in T n 
in which a is the longest flipped diagonal] for each flipped diagonal b, 
distinct from a, we have d n (b) < d n (a). 

Let T n (0) be the set of coverings in which no monomino is flipped on 
its longest diagonal. Note the distinction between a monomino flipped 
on its longest diagonal, and the longest flipped diagonal in the whole 
covering. 

The sets T n (0) and T n (a), for each diagonal a defined above, are a 
partition of T„, and the allowable diagonal flips of each subset can be 
applied independently of the other flips, by Theorem 2 in [2]. 
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4. Enumeration 

Let S(s, k) be the number of subsets of {1,2, ... ,s} whose sum is 
k. The number of coverings with k vertical (or horizontal) dominoes 
is expressible in terms of this function by making independent flips 
of diagonals whose lengths are some subset {1, 2, ... , s}. We identify 
these sets of diagonals in the proof of Lemma [TJ 

Lemma 1. Let V(n, k) and H(n, k) be the number of coverings in T n 
with exactly k vertical and horizontal dominoes, respectively. If n is 
even, then V(n, k) is equal to 



(4a) \H{n,k):=2j2 E S(n- i - 2,h)S(i - l,k 2 ) 



i=\ I ki+k 2 = 
\ ft— (n— i— 1) 



/ 



(4b) 



+ E s 

ki+k 2 =k 



n-2 



,h)S 



n-2 



When n is odd, \H(n, k) is equal to H(n, k). 



some diagonal a, and the term (4b) counts those in T n (0). 



Proof. Each outer sum term of (4a) adds the coverings for T n (a), for 



Case n even: The trivial covering in T n consists only of horizontal 
dominoes, and flipping the diagonal a contributes d n (a) vertical domi- 
noes. Diagonals l\ and r\ have even length, for all i, while I'- and r\ 
have odd length. We use this fact to find sets of diagonals which have 
lengths 1,2, ... ,s, for some s £ N, by combining allowable diagonals 
in opposite corners, for each T n (a). Table |2J shows the lengths of the 
longest allowable diagonals in each corner for each T n (a), and from this 
we can find the required sets of diagon als. F or example, the allowable 
diagonals in T n (/J) are shown in Fig. 4(a) (for (n,i) = (18,5)) and 
their respective lengths are 
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We have d n (lj) = n — i — 1, so we are interested in the number 
of combinations of the above independently flippable diagonals with 
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Figure 4. Allowable diagonals shown in alternating 
grey and white, (a), for T n (^), where (n,i) = (18,5), 
and (b), for Ti 8 (0). 



exactly k — (n — i — 1) vertical dominoes. That number is 
J2 S(n-i-2,ki)S(i-l,k 2 ). 

ki+k2 = 
k—(n—i—l) 

The indices of the diagonals Ij for which d n (lj) > d n {lf) and rf for 
which d n (rj) > d n (rj), range from 1 to [^^Jj as required for (4a). 

Now suppose a = 0. If i is the largest index such that d n (lf) < d n (lj) 
and j is the largest index such that d n {r]j) < d n (rj), then max(z, j) = 

l/^l an d K — il = 1- The allowable diagonals in T n (0) and their 
respective sizes are shown in the table below (see Fig. |4(b)). 
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Choosing subsets of the independently flippable diagonals with k 
vertical dominoes contributes the term 
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and since n - (|_(n - 2)/2j - 1) - 3 = [(n - 2)/2j , this is equal to (j4bj) 
for even n. 
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n even I- 



Index and size of largest diagonal in this corner 
(j odd) lj (j odd) rj (j even) r- (j even) 
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Table 2. The longest allowable diagonals in each of four 
corners for each T n (a). Entries are calculated using the 
parity of i and j, the avoidance of conflicts, and the re- 
quirement that a be the longest diagonal in T n (a). Note 
that "c. 2(a)", above, refers to conflict 2(a) which occurs 
between diagonals a and b if d n (a) + d n {b) > n. 



Case n odd: The trivial covering is a vertical bond with \{n — 2)/2\ 
monominoes at the top (besides the two that are fixed) and \{n — 2)/2] 
non-fixed monominoes along the bottom boundary When diagonal a 
is flipped, d n (a) horizontal dominoes are added to the covering, instead 
of vertical dominoes. Hence we argue for H(n, k) rather than V(n, k). 

Now tj~ and £+ have even length, and 6*~ and bj* have odd length 
(see Table [2]). For example, the allowable diagonals in T n (t+) are 
shown in Fig. 5(a) (for (n,i) = (17,6)), and their respective lengths 
are 
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Figure 5. Allowable diagonals shown in alternating 
grey and white, (a), for T n (tj), where (n,i) = (17,6), 
and (b), for T 17 {0). 



Once again d n (t~*) = n — i — 1, so we are interested in the number 
of combinations of the above independently flippable diagonals with 
exactly k — (n — i — 1) horizontal dominoes. As before, that number is 



y S(n — i — 2, k{)S(i — 1, &2 



k— (n— i— 1) 



Now suppose a = 0, then if z is the largest index such that d n {tj*) < 
d n (ti~) and j is the largest index such that d n {b~^) < d n (by) then 
max(i,j) = Ij^l and \i — j\ = 1. The allowable leftward diagonals in 
T n (0) and their respective sizes are given in the table below. 



i<2 , ^4 , • • • , i/j a, <±, . . . , j, 

b^ ,bf , . . . ,bf 1,3, ... ,i 

and by horizontal symmetry, the rightward diagonals have the same 
lengths. We conclude that the coverings with k horizontal dominoes of 
T n (0) is also generated by (4b) when n is odd. □ 

The terms \ffl(n, k)z k can be summed over k to obtain the generating 
polynomial T(n, z) (same as \H n (z)), mentioned in Conjecture 4 of [TJ. 



12 ERICKSON AND RUSKEY 



Theorem 1. The generating polynomial for (4a 4b) is 

(5) W n (z) : = 2 J2 S n . l ^{z)S l ^{z)z n -'- 1 + (s^j (z] 

i=l 

where S n (z) = J2k&z ^( n ' k)z k . This "generates" V(n,k) for even n, 
and H(n, k) for odd n. 

Proof. This follows from Lemma [TJ □ 

The degree of \H n (z) is IL -^ 1 — (n — 1), because this is the largest 
number of vertical dominoes possible in a covering of T n , for even n 
(and horizontal dominoes for odd n). For example, the covering with 
all U flipped up and all r t flipped down contains exactly n—1 horizontal 
dominoes. 

The coefficients of VH n (z) are listed in Table M up ton = 10, and 
the following conjecture is true at least up to n = 20. If Q(z) is a 
polynomial, then write (z k )Q(z) to denote the coefficient of z k . 
Conjecture 1. 

(a) For k < n - 2, we have (z k )W n {z) = (z k ) Umyoi 1 + zm ) 2 > the 
number of partitions of k into distinct parts with two types of each 
part (see A022567 in ©j. 

(b) For < k < n — 3, we have 

(z degiWn(z »- k )W n (z) = 2{z k ) Y[ (1 + z m ) } 

m>0 

twice the number of partitions of k into distinct parts (see A000009 
in \§\). 

Rotating a covering of T n by 7r/2 radians interchanges vertical and 
horizontal dominoes, and this transformation can be applied to the 

2_ 

generating polynomial \H n (z) to obtain the polynomial \H n (z~ 1 )z^^ . 
Thus we can easily derive the bivariate generating polynomial R n {x, y), 
whose coefficient of x v y h is the number of tatami coverings with exactly 
v vertical dominoes and h horizontal dominoes. 
Our remarks prove the following corollary. 

Corollary 1. Let R n (x,y) be as defined above. We have 

(6) Rn(x,y) = 2W n (xy- 1 )y rJ ^ 1 + 2W n (x- 1 y)x^ . 

We list some basic properties of R n (x,y). 

• The degree of R n (x, 1) as well as the degree of every term in 
R n {x,y) is ^f^; 
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n\z k 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 



2 


1 
































3 


1 


2 






























4 


1 


2 


3 


2 


























5 


1 


2 


3 


6 


4 


2 


2 




















6 


1 


2 


3 


6 


9 


8 


7 


6 


2 


2 


2 












7 


1 


2 


3 


6 


9 


14 


15 


14 


11 


10 


8 


6 


1 


2 


2 


2 


8 


1 


2 


3 


6 


9 


14 


22 


24 


25 


28 


25 


22 


19 


14 


10 


10 


9 


1 


2 


3 


6 


9 


14 


22 


32 


37 


42 


49 


18 


49 


46 


38 


31 


10 


1 


2 


3 


6 


9 


14 


22 


32 


16 


56 


66 


78 


81 


90 


92 


88 


n\z k 


16 


17 


18 


19 


20 


21 


22 


23 


24 


25 


26 


27 


28 


29 


30 


31 



8 


8 


4 


4 


2 


2 


2 






















9 


30 


21 


20 


16 


12 


12 


10 


6 


4 


4 


2 


2 


2 








10 


81 


76 


69 


58 


51 


11 


38 


31 


28 


22 


20 


16 


14 


12 


8 


6 


n\z k 


32 


33 


34 


35 


36 


37 


38 


39 


40 


41 


42 


43 


44 


45 


46 


47 



10 4 4 2 2 2 

Table 3. Table of coefficients of \H n {z) for 2 < n < fO. 
The (n, /c)th entry represents the number of coverings 
of T n with k vertical dominoes when n is even, and k 
horizontal dominoes when n is odd. 



• the polynomial R n (x,y) can be recovered from R n (x, 1), and 
the latter is the generating polynomial for the set of all n x n 
coverings with n monominoes with exactly v vertical dominoes 
(or h horizontal dominoes); 

• the polynomial R n {x, 1) is self reciprocal because of interchange- 
ability of vertical and horizontal dominoes; and finally, 

• the polynomial R n (x,l) has similar properties to those listed 
for \H n (z) in Conjecture [II in the sense that for some increas- 
ing integer function /, we have (x k )R n (x,l) = (x k )R n+1 (x, 1), 
whenever k < f(n + 1). 

If there is an even number of dominoes, which is the case when 
{n 2 — n)/4 is an integer, then (x k y k )R n (x,y) = 4(z k )\H n (z) , where 
k = (n 2 — n)/4. Rotating the covering maps k vertical dominoes to k 
horizontal dominoes, and vice versa. The coverings counted by these 
coefficients are called balanced tatami coverings, appropriately named 
by Knuth (private communication), because the number of vertical and 
horizontal dominoes are equal. Here is (z k )\tl n (z) for 2 < n < 56: 0, 
0, 2, 2, 0, 0, 10, 20, 0, 0, 114, 210, 0, 0, 1322, 2460, 0, 0, 16428, 31122, 
0, 0, 214660, 410378, 0, 0, 2897424, 5575682, 0, 0, 40046134, 77445152, 
0, 0, 563527294, 1093987598, 0, 0, 8042361426, 15660579168, 0, 0, 
116083167058, 226608224226, 0, 0, 1691193906828, 3308255447206, 
0, 0, 24830916046462, 48658330768786, 0, 0, 366990100477712, (see 
A182107 in |6J). Note that this is perhaps better viewed as four se- 
quences, one for each < j < 4 such that n (mod 4) = j. 
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5. Exhaustive generation of tatami coverings in constant 

amortized time 

In this section, we present procedure genVH(n, k), which generates 
the coverings counted by \H(n,k), while doing a constant amount of 
data structure change per covering that is produced (example output is 
shown in Fig. [6]). Let S(n, k) denote the set of subsets of {1, 2, ... , n} 
whose elements sum to k\ thus |S(n, k)\ = S(n,k). The procedure 



follows naturally from the sums in Equation (4a 4b), since each term 
S(a,i)S(b,j) counts some set of '/' and '\'-oriented diagonals. The 
sets S(o, i) x S(b,j) are generated in constant amortized time (CAT) by 
a modification of C4 from [7] (see Listing [I]). Our modified algorithm, 
modC, is invoked for each sum term of Equation (4a 4b). Procedure 
modC is CAT for the same reasons that C4 is CAT. 
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Figure 6. The coverings of T§ with exactly 7 vertical 
dominoes. This is the output of genVH(8, 7) printed in 
the order the coverings are generated (as one would nat- 
urally read text). 



There is one subtlety involved in exploiting the CATness of C4. In- 
voking C4(a, i) requires Q(a) preprocessing steps if its input list is 
recreated for each call, but C4(a, i) may not produce so many com- 
binations for small a and large i. The result is that we may make 
many calls to modC that require too much preprocessing, but this is 
dealt with, as follows: a top level call to C4(i,j) in [7] takes the list 
[i+l, 2, 3, ... , i+1], which requires i+1 steps to create, however, CA(i,j) 
also concludes with the same list (see Listing [I]) . Let A and B be the 
largest integers for which modC is called to compute S(a,i) x S(b,i). 
We set aiSet= [1, 2, 3, . . . , A] and biSet= [1, 2, 3, . . . , B], and by set- 
ting aiSetfO] = a + 1 and biSet[0] = b + 1, we initialize for each call 
to modC with exactly two operations. 
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Listing 1. Python code for a modified version of 
C4 from j7] to compute S(o, i) x S(b,j). Global vari- 
ables aiSet and bjSet are the lists representing S(a, i) 
and S(6, i), respectively 
def modC (a , i , b , j , comp , isFirst ) : 
global aiSet,bjSet 
if ( a == 0) : 

if ( isFirst ) : 

modC(b,j ,0,0, False , False) 
else : 

Output (aiSet , b j Set ) 
else : 

if ( isFirst ) : 

L = aiSet 
else : 

L = bjSet 
if ( i > a*(a+l)/2 ) : 

i = a*(a+l)/2 - i; comp = not comp 
if ( i<a ) : 

if ( comp ) : 

L[a] = L [0] ; L [0] = i + 1 
modC ( i, i, b, j, comp, isFirst) 
L [0] = L[a] ; L [a] = a+1 
else : 

modC ( i, i, b, j, comp, isFirst) 
else : 

L[a] = L [0] ; L [0] = a 
if ( comp ) : 

modC ( a-1 , i, b, j, comp, isFirst) 
L [0] = L[a] ; L [a] = a + 1 
modC ( a-1, i-a, b, j, comp, isFirst) 
else : 

modC ( a-1, i-a, b, j, comp, isFirst) 

L [0] = L[a] ; L [a] = a + 1 

modC ( a-1, i, b, j, comp, isFirst) 



Theorem 2. The coverings in T n with exactly k vertical dominoes 
if n is even and horizontal dominoes if n is odd, can be exhaustively 
generated in constant amortized time. 
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Proof. The outer procedure does a constant amount of work per call to 
modC. This subroutine is CAT, so the outer procedure is also CAT. □ 



6. A MYSTERIOUS FACTOR OF W n (z) 

In this section we prove that the generating polynomial \H n (z) has 
(very nearly) the factorisation conjectured in pQ. We use the following 
lemma. 



Lemma 2. For all x > 0, 

(7) [*J = E 



k>l 



X 1 

¥ + 2 



Proof. Let n = \_x\ and apply strong induction on n. Clearly Equa- 
tion ([7| holds for the base case, when n = 0. Suppose it holds for 
0, 1 Ti — 1, then 



fc>i 



X 1 

¥ + 2_ 


= 


x 1 

_2 + 2_ 


fc>i 


2 
2 fc 


1 

+ 2. 










= 


x 1 

_2 + 2_ 


-LIJ 


= 


2 


1 
+ 2. 


+ 


_ 2 



Lxj. 



Two applications of Equation (3.11) in [3] yield the penultimate equa- 
tion and the final equation follows by considering the parity of [x\ , or 
by using (3.26) in [3] with m = 2 and x' = x/2. D 

Theorem 3. The generating polynomial VH n (z) has the factorisation 

W n (z) = P n (z)D n (z) 

where P n {z) is a polynomial and 



DJz) 



Us 

i>i 



L^J 



[z). 



Proof. We prove that D n (z) divides \H n (z) by using the factorisation 
of S n (z) into cyclotomic polynomials ([TJ, Lemma 5), 



(9) 



S n (z) = H^(z) [ 



n+1 
2j - 



and showing that the power of $i(z) is greater in each term of VH n (z) 
than it is in D n (z). 
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The power of &2j(z) in D n (z) is obtained by substituting Equa- 
tion ([9]) into Equation (|8j, 



2/ 



i>l 



»>1J>1 



z2i>l 



2j 



+J 



We simplify D n (z) to 
(10) !>„(«) =J]*^) L ^- 1 , 

by applying Lemma |2J and with Equation (3.11) in [3]. 
Expanding the second term of \H n (z) gives 



which is divisible by D n (z), since 

n-2 

< 2 



-+j 



2.; 



n-2 



2j 



+ J 



2j 



for all j > 1 and positive even integers n. 
The other terms in V n (z) are of the form 



(n-fc-2)+j 
2j 



j) (nw*^) 



5 n _jfc_ 2 (^5 fc _i^)^ = n^*' 

\i>o 

for each 1 < fc < [^y^J where d is the appropriate power of z. These 
terms are all divisible by D n (z) if the exponents in Equation ( 10 ) satisfy 



(11) 



n-2 



< 



k-1 1 
~^T + 2 



+ 



n-k-2 1 

2j + 2 



Let ri and r 2 be integers such that < r^ < 2j and 



fc-i 

2j 



fc-i 



and ^2 

2j 



n-2 
2j 



+ ^. We eliminate occurrences of 



from Inequality (11), since they are integers and can 



2.7 
and 



fc-i 
2i 

be removed 



^ 2j 
n-2 
2 -? . 

from 



floors, and rewrite the inequality as 



(12) 



0< 



T\ 1 




77T + " 


+ 


L 2 J 2 . 





r 2 — ri — 1 1 

2j + 2 
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It is straightforward to show that if the second term is —1, then the 
first term is equal to 1. 

Therefore, D n (z) divides each and every term of VH n (z). □ 



n\z 
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10 


11 


12 


13 


14 


15 
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2 
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1 


1 


2 
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1 
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2 


4 





2 






















6 


1 
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2 
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-2 


2 




















7 


1 





1 


2 
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4 


-2 


4 





2 


-2 


2 










8 


1 





1 


1 


2 


3 


1 


-2 


2 





4 


-2 


2 


-2 


2 




9 


1 





1 


1 


2 


3 


1 


(5 


-2 


6 





8 


-2 


1 


-1 


6 


10 


1 


-1 


1 





1 


1 


1 


2 


2 


-6 


6 


-2 


6 


-6 


4 


-1 


11 


1 


-1 


1 





1 


1 


1 


2 


2 


1 


-8 


10 


-4 


10 


-8 


8 


n\z k 


16 


17 


18 


19 


20 


21 


22 


23 


24 


25 


26 


27 


28 


29 


30 


31 


9 


-2 


4 


-2 


2 


-2 


2 






















10 


6 


-6 


6 


-4 


4 


-4 


2 




















11 


-8 


10 


-10 


12 


-8 


10 


-12 


10 


-6 


6 


-6 


6 


-1 


1 


-1 


2 



Table 4. Table of coefficients of P n (z) for 3 < n < 11. 



Our computer investigations show that P n (^) is irreducible for 1 < n < 200, 
and we know the complete factorisation of Sk(z), for each positive in- 
teger k. We suspect, therefore, that the complete factorisation is 

(13) W n (z) = P n (z) H $ 2j (z) L^J . 

i>i 

The factor P n (z) is somewhat more mysterious than D n (z); e.g., we 
have no formula to express it besides \H n (z)/D n (z). Take P\\{z) for 
example, which is equal to 1 — lz 1 + lz 2 + Oz 3 + lz 4 + lz 5 + lz 6 + 
2^ 7 + 2z 8 + Az 9 - 8z w + lOz 11 - Az 12 + 10z 13 - 8z 14 + 8z 15 - 8z w + 
10z 17 - 10z 18 + 12z 19 - 8z 20 + 10z 21 - Viz 22 + IO2 23 - 6z 24 + 6z 25 - 
Qz 26 + 6z 27 - Az 28 + Az 29 - Az 30 + 2z 31 . The coefficients are almost all 
non-zero, the central coefficients are larger than the ones at the tails, a 
great many of them are even, they alternate in sign for a long stretch, 
and the polynomial is irreducible. 

The degree of P n (z) is deg(P n (^)) = deg(W n (z)) - deg(D n (z)), 
both of which are easily calculated. In general deg(P n (z)) is equal to the 
sum of the sequence of largest odd divisors of the numbers 1, 2, . . . , n — 
2, which is a sequence with some nice properties (see A135013 in [6]). 

Theorem 4. For each n > 2, 

n-2 

deg(P n (z)) = J20d(k), 

k=l 

where Od(k) is the largest odd divisor of k. 
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(14) 



deg(P n (z)) 



since deg(S n (z)) = ( n J 1 ). 

The proof that Ylk=i Od(k) = deg(P n+ 2(-2)) is by induction, and the 
base case, where n = 0, is easily verified. Let p n = deg(P n (,z)), for 
n > 2, to abbreviate the notation. It remains for us to show that 
Pn+z-Pn+2 = Od{n + l). 

Let n'01 a be the binary representation of n, so that n + 1 = n'10 a , 
and let {AJ = 1 if the statement A is true, and 
Observe that 



otherwise. 



(15) 



71 + 1 



2 A ' 



[*<«] + 



n 

2k 



which we use to simplify 



m [ n +i 



k>l 



L^J + i 

2 



and write 



a 

Pn+3 ~ Pn+2 = (« + 1) - ^ ( [ok 

k=l 



+ 1 • 



Using Equation (15) and the fact that (n + l)/2 is an integer for 
1 < k < a, we write 

and then express this as the remaining sum terms in Equation (J7l) 

Pn+3 - Pn+2 



, A / n + 1 1 \ 

k=1 \ L J / 



E 

fc>ct+l 

E 

fc-a>l 



n + 1 1 

2 k + 2 



71+1 

2 a 



2 fc-a + 2 



Applying Equation ([7]) again, we have p n+ 3 — p n +2 = (n + l)/2 a , 
which is equal to Od{n + 1), as required. □ 

In addition to finding deg(Pn(;z)), we can evaluate at z = 1 with 
Pn(l) = Wn(l)/Pn(l) = 22, a ratio which is also easy to calculate 
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in general because W n (l) and D n {l) have well understood combina- 
torial interpretations. It also leads to an interesting sequence, whose 
derivation for all n is given below. 

Theorem 5. The sum of the coefficients ofP n (z) is equal to n2 v ^ n ~ 2 ^~ x , 
where u(n) is the number of Is in the binary representation of n 

Proof. The sum of the coefficients of P n (z) is equal to P n (l), which is 
expressible as \H n (l) / D n (l) . The numerator evaluates to n2 n ~ 3 , since 
this is the number of coverings in T n , and the denominator is evaluated 
as described below. 

It is well known that $fe(l) = p if k is a non-zero power of a prime 
p and $fc(l) = 1 if k is divisible by two distinct primes (see jl], p. 74). 



We can evaluate D n (l) using Equation (10), 



D n {\) = Jl $ 2i (l) L^J = 2^ L^J , 



»>1 

by ignoring the factors for which 2% is not a power of 2. Apply Equation 
(4.24) in to obtain D n (l) = 2 n - 2 -^ n ~ 2 \ Thus 

P n (l) = n 2 n - 3 - (n " 2)+i/( "" 2) = n2 u(n - 2) -\ 

D 
We have verified that P n {z) is irreducible over the integers for 1 < n < 200, 
but we do not understand its structure well enough to prove it for all n. 
We state below some of the observable structure which has also been 
verified for 1 < n < 200, as Conjecture [2j and we plot some complex 
roots for odd n up to 67 in Fig. [7| 

Conjecture 2. 

(a) If k > 1 and n (mod 2 k ) = 2, then (z^P^z) = (^)P n+j (^) for 
i<§^ andj<2 k . 

(b) When n is odd, P n (z) has exactly one real root a n , with — 1 < a n < —0.5, 
and {a n } n m is a monotonically decreasing sequence. 

(c) When n is even, P n (z) has no real root. 

(d) The polynomial P n (z) is irreducible over the integers for n > 2. 

(e) The alternating sums of coefficients are given by the generating 
function 



(16) 



yp.(-i).--'= ( l + ^ l zli 



n>2 



1 -2z 2 )Vl ~4:Z 2 ' 



(f) For even n, the sum of the absolute values of coefficients of P n (z) 
is equal to P n (— 1) when n > 20. 
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Figure 7. The complex zeros of P n (z) for odd 
n, where 3 < n < 67. Darker and smaller 
points are used for larger n. Larger versions may 



be viewed at http://webhome.cs.uvic.ca/~ruskey/ 
Publications/Tatami/HoriVert . html. 



The right hand side of Equation (16) is the sum of two generating 
functions, with odd and even powered terms, respectively. The se- 
quence of coefficients of the odd power terms is — ^2 i=0 2 k ~ l ( /) , for 
k > (see A082590 in [6]), and that of the even power terms is („), 
for k > 1 (see A000984 in [6]). The first few numbers P n {— 1), starting 
with n = 2, are: 1, -1, 2, -4, 6, -14, 20, -48, 70, -166, 252, -584, 
924, -2092, 3432, -7616, 12870, -28102, 48620, -104824, 184756, 



-394404, 705432, 



-1494240, 2704156, 



-5692636, 10400600. 
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Conjecture |2j[fJ) compares the above sequence with the sum of the 
absolute values of the coefficients of P n (z). The first few of these are 
listed, also starting with n = 2: 1, 3, 4, 10, 10, 22, 28, 64, 76, 180, 
260, 606, 932, 2124, 3440, 7666, 12872, 28178, 48620, 104946, 184756, 
394638, 705432, 1494600, 2704156, 5693376, 10400600. 

7. Conclusions and further research 

The polynomials P n {z) exhibit numerous patterns in the signs of 
their coefficients, their plots and zeros (see http : //webhome . cs . uvic . 
ca/~ruskey/Publications/Tatami/HoriVert .html), their degrees, and 
the values of P n (—1), but yet all we have seen is that they fall mag- 
ically out of these coverings. What is the geometric interpretation, if 
there is one, for the factorisation of \H n (z), and how do we calculate 
the coefficients of P n (z) without dividing D n (z) into \H n {z)l 

In the present paper we deal with nxn coverings with n monominoes, 
and the techniques can be readily applied to r x c coverings with r < c 
and a maximum number of monominoes. When m is not maximum, 
however, we may have to deal with other features of tatami coverings, 
called bidimers and vortices, described in [TJ. The numbers of hor- 
izontal and vertical dominoes that any given bidimer or vortex will 
introduce is easily calculated, and they provide the advantage of some- 
times isolating corners of the grid, which makes diagonal flips easier to 
count. We may not, however, have the good fortune of encountering 
a formula like the one in Lemma [TJ which easily yields the generating 
polynomial of Theorem [TJ 

A generating function for fixed height balanced coverings, or perhaps 
some other relation between the numbers of vertical and horizontal 
tiles, would extend the results on fixed height coverings in [3] and pQ. 

Acknowledgements: We thank Prof. Don Knuth for his valuable 
communications about tatami coverings. 
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